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Abstract 


Adamantanes (amantadine and rimantadine) have been used to prevent and treat influenza 
A virus infections for many years; however, resistance to these drugs has been widely re- 
ported in the world. To investigate the frequency and distribution of M2 gene mutations in 
adamantane-resistant influenza variants circulated in the world between 1902 and 2013, 
31251 available M2 protein sequences from different HA-subtype influenza A viruses (H1— 
H17) were analyzed and adamantane resistance-associated mutations were compared 
(L26F, V27A, A30T, A30V, S31N, G34E, and L38F). We find that 45.2% (n = 14132) of influ- 
enza A (H1—H17) viruses circulating globally were resistant to adamantanes, and the vast 
majority of resistant viruses (95%) bear S31N mutations. Whereas, only about 1% have 
V27A mutations and other mutations (L26F, A30T, G34E, and L38F) were extremely rare 
(their prevalence appeared to be < 0.2%). Our results confirm that H1, H3, H5, H7, H9, and 
H17 subtype influenza A viruses exhibit high-level resistance to adamantanes. In contrast, 
the appearance of adamantane-resistant mutants in H2, H4, H6, H10, and H11 subtypes 
was rare. However, no adamantane resistance viruses were identified among other HA sub- 
types (H8, H12—H16). Our findings indicate that the frequency and distribution of adaman- 
tane-resistant influenza variants varied among different HA subtypes, host species, years 
of isolation, and geographical areas. This comprehensive study raises concerns about the 
increasing prevalence of adamantane-resistant influenza A viruses and highlights the im- 
portance of monitoring the emergence and worldwide spread of adamantane- 

resistant variants. 


Introduction 


Influenza A virus is a common cause of respiratory infections, infecting humans, birds, pigs, 
horses, as well as other species [1]. Eighteen Hemagglutinin (HA, H1 to H18) and eleven 
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Neuraminidase (NA, N1 to N11) subtypes of influenza viruses can cause epidemics and pan- 
demics through antigenic drift and antigenic shift, respectively [1]. Considering the significant 
threat that influenza A viruses pose to global economy and human health, preparedness for a 
potential influenza pandemic is a global priority [2]. 

The primary means to protect humans against influenza virus infection is by vaccination. 
However, in the absence of an effective and widely available vaccine, protection from a pan- 
demic and/or long-term endemicity would rely largely on the prophylactic and therapeutic 
properties of antiviral treatment. Therefore, effective anti-influenza drugs offer a valuable addi- 
tion to the use of vaccines, particularly for those not caused by drug-resistant viruses [3]. Cur- 
rently, two classes of anti-influenza drugs, adamantanes (amantadine and its methyl derivative 
rimantadine) and neuraminidase inhibitors (oseltamivir and zanamivir), are used for prophy- 
laxis and treatment of influenza A infections [4]. These FDA-approved drugs are effective for 
influenza A viruses, including highly pathogenic H5N1 viruses [5]. Their therapeutic efficacy, 
however, has come into question with the emergence of drug-resistant viruses. Adamantane 
drugs, also known as M2 channel blockers, could block the ion channel formed by the M2 pro- 
tein encoded by the M gene of influenza A viruses and inhibit early stages of virus replication 
[6]. As adamantane derivatives, amantadine (SymmetrelTM) and rimantadine (FlumadineTM) 
have been used successfully for influenza A virus infection for more than 30 years because of 
their wide availability and low cost [7, 8]. The effectiveness of prophylaxis with these drugs var- 
ied between 80% and 90%, and their use can reduce the duration of illness by 1.5 days if treat- 
ment is begun within 48 hours of infection [9]. However, the main drawbacks of ultilizing 
adamantane derivatives are that drug-resistant variants develop rapidly [10]. Presently, animal 
experiments and clinical studies have indicated the frequent occurrence of amantadine-resis- 
tant influenza variants after exposure to the drugs. During clinical use in humans, adamantane 
resistance emerged in about 30% of patients after 2-3 days of the treatment [11]. Moreover, 
drug-resistant variants can replicate as efficiently as sensitive ones, and they can also transmit 
efficiently from one individual to other individuals without apparent loss of pathogenicity [11]. 
Additionally, complete cross-resistance exists between amantadine and rimantadine because of 
the common mechanism of action [12]. To reduce the emergence and spread of drug-resistant 
influenza viruses, amantadine-oseltamivir combination chemotherapy has been evaluated and 
their synergistic antiviral effects have been confirmed [13]. Therefore, understanding the 
mechanism responsible for the emergence and spread of antiviral resistance is important for 
controlling seasonal and pandemic influenza. 

After four decades of effective use of adamantanes in the prophylaxis and treatment of 
influenza, global resistance to these drugs has increased dramatically. In the period from 1991 
to 1995, the global frequency of adamantane resistance among A (H3N2) viruses was as low as 
0.8% [14]. However, in recent years, a significant worldwide increase in adamantane resistance 
has been reported for human H1N1 and H3N2 seasonal influenza viruses and H5N1 avian in- 
fluenza viruses [15-17]. For example in China, the frequency of amantadine-resistant HIN1 
variants has greatly increased from 28% during 2004-2005 to 72% in 2005-2006 [18]. Similar- 
ly, in Taiwan and Canada, the frequencies of adamantane-resistance raised from 0 in 2004- 
2005 to 24% and 33% in 2005-2006, respectively [18]. In the United States, resistance to ada- 
mantine among A (H3N2) viruses circulating had reached 92% in 2005 [19]. Specially, ada- 
mantane resistance had reached 100% during 2005-2006 in influenza A (H3N2) isolates from 
some Southeast Asian countries [18]. It is worth mentioned that the most recent H5N1 strains 
isolated from humans and poultry in Vietnam and Thailand show resistance to the adaman- 
tanes [20]. The rapid emergence and spread of adamantane-resistant variants of influenza 
H1N1, H3N2, recent H7N9 viruses from China, and H5N1 viruses have led to changes in CDC 
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recommendations for the use of adamantanes in the control of influenza A virus infections 
[21]. 

The genetic basis for resistance to adamantanes has been well characterised and amino acid 
substitutions at residues L26, V27, A30, S31, and G34 within the transmembrane domain of 
the M2 protein are associated with adamantane resistance [12, 22]. Abed et al. constructed re- 
combinant H1N1 influenza A viruses with the commonly observed L26F, V27A, A30T, S31N, 
G34E, and V27A/S31N mutations in the M2 gene by using reverse genetics, and showed all of 
these mutations are able to cause amantadine resistance [23]. So far, the known mutations that 
confer resistance against adamantanes are L26F, V27A, A30T, A30V, S31N, G34E and L38F in 
the M2 transmembrane region [6, 22]. Bright et al. found that out of 7000 influenza A isolates 
collected from countries in Africa, Asia, Europe, the Americas and Oceania between 1994 and 
2005, the frequency of adamantane resistance mutations rose from 0.4% in 1994-1995 to 
12.3% in 2003-2004 [16]. Most drug-resistant influenza viruses contain one of these amino 
acid changes, but variants with dual mutations have also been described [12]. In particular, the 
S31N substitution renders the virus highly resistant to adamantanes, and the incidence of clini- 
cal isolates bearing this mutation has jumped from a small percentage to approximately 97% in 
recent years [18, 19, 24]. Given the uncertainty of the therapeutic value about adamantanes, 
particular to a pandemic and/or long-term endemicity, it is necessary to study deeply whether 
resistance to adamantanes is prevalent in various influenza viruses in the world. In this present 
study, we analyzed the incidence rate, host distribution, temporal distribution, and geographic 
distribution of adamantane-resistance among H1-H17 subtype influenza A viruses isolated 
worldwide from 1902 to 2013, and also tried to clarify the possible reasons of high adaman- 
tane-resistance incidence. 


Methods 
Sequence data 


All amino acid sequence data of the M2 proteins used in this study were obtained from the in- 
fluenza virus sequence database available in GenBank (Influenza Virus Resource [http://www. 
ncebi.nlm.nih.gov/genomes/FLU/FLU. html], accessed on 8 June and 31 October 2013). Firstly, 
the shorter and identical sequences were collapsed using the web servers of Influenza Virus Re- 
source in the National Center for Biotechnology Information (NCBI), and the remaining se- 
quences were then utilized for further analysis. Secondly, a total of 31251 M2 protein 
sequences (nucleotide positions 1-27 and 716-979) were selected and downloaded represent- 
ing the spectrum of influenza A diversity from 1902 to 2013, together with information about 
the subtype, host, location and year of isolation. Subtypes H1-H17 of influenza A virus from 
different host species, including human, avian, swine, equine, and canine viruses, were includ- 
ed. Asa result, M2 protein sequences of 12579 H1, 539 H2, 9414 H3, 1001 H4, 2856 H5, 1208 
H6, 1147 H7, 124 H8, 1151 H9, 548 H10, 421 H11, 133 H12, 74 H13, 11 H14, 7 H15, 35 H16, 
and 3 H17 influenza A viruses were used for the analysis. All GenBank accession numbers of 
adamantane resistant viruses were listed in the S1—S6 Tables through primary classification in 
host, subtype, and isolation location. 


Sequence Analysis 


All of these sequence data were aligned by Clustal W (BioEdit version 7.0.5) (http://www.mbio. 
ncsu.edu/BioEdit/bioedit.html). Substitutions of residues Leu26, Val27, Ala30, Ser31, Gly34, 
His37, Leu38, and Trp41 in the transmembrane region of the M2 ion channel protein were 
used to screen amantadine-resistant mutants [12, 25]. The frequencies of the appearance of 
recognized adamantine resistance-associated amino acids (L26F, V27A, A30T, A30V, S31N, 
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G34E, and L38F) were identified and tabulated. Among 31251 sequences, the number of strains 
with mutations was shown, and their detail information was listed in the S1-S6 Tables. 


Results 
The frequency of adamantane resistant influenza variants 


A total of 31251 influenza A viruses were screened for the specific mutations known to corre- 
late with resistance to adamantanes. Our analysis based on the M2 protein sequence data avail- 
able in the GenBank showed a diversity in the frequency of adamantane resistant influenza 
variants among different HA subtypes, host species, years of isolation, and geographical areas 
(S1-S6 Tables). In the present study, a total of 14132 (45.22%) of the 31251 influenza A viruses 
contained one or two amino acid substitutions in the M2 protein known to cause resistance to 
adamantanes [12, 22]. The prevalence of adamantane resistance over the time period of the 
study is shown in Table 1. These results suggested that the frequency of emergence of drug-re- 
sistant strains varied among different HA subtypes. Among these isolates, we identified ada- 
mantane-resistant influenza variants from H1-H7, H9-H11, and H17 subtypes. However, 
analysis of viruses isolated from 1902 to 2013 did not reveal H8 and H12-H16 subtype influen- 
za isolates with amino acid substitutions in the transmembrane region of M2 protein corre- 
sponding to resistance to adamantanes (Table 1). In other words, no adamantane resistant 
viruses were identified among the influenza H8, H12, H13, H14, H15, and H16 subtype. In 
contrast, adamantane-resistant strains frequently occurred among viruses of six HA subtypes 
with the frequencies of H1 (8777/12579, 69.78%), H3 (4107/9414, 43.63%), H5 (804/2856, 
28.15%), H7 (144/1147, 12.55%), H9 (265/1151, 23.02%), and H17 (3/3, 100%). The remaining 
HA subtype influenza viruses were identified with <0.5% adamantane resistant variants except 
H6 subtype influenza viruses. Analysis of sequence data from H6 subtype strains indicated 
2.07% (25/1208) resistant variants with M2 protein-associated mutations that would confer re- 
sistance to adamantanes [22]. Notably, the numbers of resistant variants differed among HA/ 
NA subtype (Table 1). For example, a total of 8193 influenza A (H1N1) viruses and 583 influ- 
enza A (H1N2) isolates were found to be adamantane-resistant. Similarly, 4095 influenza A 
(H3N2), 7 influenza A (H3N1) and 5 influenza A (H3N8) viruses were found to contain muta- 
tions conferring resistance to adamantanes (Table 1). Our study further suggested that the fre- 
quency of emergence of adamantane-resistant strains varied among different HA or HA/ 

NA subtypes. 

The frequencies of adamantane resistance mutations of viruses from influenza H1-H17 
subtypes are also analyzed in this study. The critical residues conferring resistance to adaman- 
tanes were shown in Table 2. Results from our analyses indicated that the viruses with L26F, 
V27A, A30T, G34E, and L38F mutations were found sporadically. However, the majority of 
drug-resistant influenza variants contained point mutations resulting in a serine to-asparagine 
change at amino acid 31 (S31N) of the M2 protein that conferred adamantane resistance [22]. 
Among these resistant viruses, 96.7% (8490/8777) of the H1, 93.6% (3844/4107) of the H3, 
66.7% (2/3) of the H4, 83.3% (670/804) of the H5, 92% (23/25) of the H6, 86.1% (124/144) of 
the H7, 87.5% (232/265) of the H9, 100% (1/1) of the H10, 50% (1/2) of the H11, and 100% 
(3/3) of the H17 strains demonstrated the S31N substitution (Table 2). Furthermore, influenza 
H1, H3 and H5 variants contained the L26F, V27A, A30T, L26F/S31N, and V27A/S31N 
adamantane resistant mutations (Table 2). These resistant strains possessed amino acid 
substitutions in the M2 protein at three positions-L26F, V27A, and S31N [22]. The G34E and 
G34E/S31N mutations emerged only in influenza A H1 viruses, whereas the L38F mutation 
only occurred in the influenza A H2 virus. In addition, the L38F/S31N and A30V/S31N 
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Table 1. The frequencies of adamantane-resistant variants among all influenza A (H1-H17) viruses used in this study. 


Viruses The number of different HA/NA subtype Percentage ofresistant 
resistant variants variants(%)° 
HA From The number of resistant Year of 
subtypes GenBank variants ° isolation? 
H1 12579 8777 1933-2013 H1N1(8193), H1N2(583) H1(1) 69.78 
H2 539 1 1968 H2N2(1) 0.19 
H3 9414 4107 1968-2013 H3N2(4095), H3N1(7), H3N8(5) 43.63 
H4 1001 3 2003, 2011 H4N2(2), H4N6(1) 0.30 
H5 2856 804 1959-2013 H5N1(766), HS5N2(26), HSN5(5), HSN3(3), 28.15 
H5(4) 
H6 1208 25 1999-2013 H6N1(20), H6N2(3), HEN6(1), H6N8(1) 2.07 
H7 1147 144 1902-2013 H7N9(65), H7N2(45), H7N7(31), H7N3(2), 12.55 
H7N1(1) 
H8 124 0 0 
H9 1151 265 1998-2013 H9N2(262), H9N8(1), H9(2) 23.02 
H10 548 1 2012 H10N7(1) 0.18 
H11 421 2 2007, 2011 H11N9(1), H11N3(1) 0.48 
H12 133 0 0 
H13 74 0 0 
H14 11 0 0 
H15 m 0 0 
H16 35 0 0 
H17 3 3 2009, 2010 H17N10(3) 100 
Overall 31251 14132 1902-2013 H1-H17(14132) 45.22 


a adamantane resistance-associated mutations present in the M2 Protein of viruses were at positions L26F, V27A, A30T, A30V, S31N, G34E, and L38F 
that confer resistance to adamantanes. 

° Appearance time to adamantane resistant variants. 

° The frequency of adamantane resistant variants among different HA subtypes. 


doi:10.1371/journal.pone.0119115.t001 


Table 2. The frequency of adamantane resistance mutations of viruses from different influenza HA subtypes. 


Adamantane-resistant mutation The number of adamantane resistance mutations of viruses from different influenza HA subtypes 
H1 H2 H3 H4 H5 H6 H7 H9 H10 H11 H17 

L26F 2 8 4 2 

V27A 5 9 95 3 26 1 

A30T 2 2 4 1 6 

S31N 8490 3844 2 670 23 124 232 1 1 3 

G34E 2 

L38F 1 

L26F /S31N 23 16 1 10 1 

V27A/S31N 250 224 21 1 15 

A30T/S31N 1 1 

A30V/S31N 2 

G34E/S31N 2 

L38F/S31N 1 

Overall 8777 1 4107 3 804 25 144 265 1 2 3 


doi:10.1371/journal.pone.0119115.t002 


PLOS ONE | DOI:10.1371/journal.pone.0119115 March 13, 2015 5/20 


IO PLOS | ONE 


Frequency and Distribution of M2 Gene Mutations 


Table 3. Host distribution of adamantane resistant influenza variants analyzed in this study among the different HA subtypes. 


The number of adamantane resistant influenza variants among different hosts overall 


Human Avian Swine Wild boar Canine Cheetah Ferret Mink Cat equine bat environment unknown 


HA subtypes 

H1 7232 
H2 1 

H3 3568 
H4 

H5 143 
H6 

H7 10 
H9 

H10 

H11 

H17 

Overall 10954 


6 


1032 


doi:10.1371/journal.pone.0119115.t003 


1525 


525 


2062 


1 


2 2 4 1 4 8777 


Hat 2 144 


mutations were only found in influenza A H3 viruses. Our study further revealed adamantane 
resistance mutations appeared across different HA subtypes. 


Host distribution of adamantane resistant influenza variants 


Adamantane resistance is caused by a single point mutation at the following amino acid posi- 
tions in the transmembrane region of the M2 protein: 26, 27, 30, 31, 34, and 38. The known mu- 
tations that confer adamantane resistance are L26F, V27A, A30T, A30V, S31N, G34E and L38F 
[12, 22]. To better understand host distribution of adamantane resistance, we conducted the M2 
protein sequence analyses focusing on positions 26, 27, 30, 31, 34, and 38 for amino acid substi- 
tutions associated with naturally occurring adamantane resistance. The numbers of adamantane 
resistant influenza variants with drug resistance mutation(s) among different hosts are shown 

in Table 3 and Table 4. Of those, adamantane resistant influenza A H1 and H3 variants mainly 
emerged in human and swine. In contrast, influenza A H5, H7 and H9 subtype viruses with ada- 
mantane resistance primarily occurred in avian. Particularly, influenza A H4, H6, H10, and H11 
subtype adamantane resistant strains were only detected in avian species, whereas influenza H2 
and H17 subtype adamantane resistant variants were only found in human and bat, respectively 
(Table 3). Adamantane-resistant mutations known as L26F, V27A, A30T, S31N, G34E, and 
L38F were detected at all six sites (positions 26, 27, 30, 31, 34, and 38) of M2 protein in human 
influenza viruses (Table 4). We also found influenza strains with double adamantane resistance 
mutations (L26F/S31N, V27A/S31N, A30T/S31N, A30V/S31N, G34E/S31N, and L38F/S31N) 

at positions 26 and 31, 27 and 31, 30 and 31, 31 and 34 or 31 and 38 in the human population. 
Adamantane-resistant mutations (L26F, V27A, A30T, S31N) were found at positions 26, 27, 30, 
and 31 in avian and swine influenza viruses. Double resistance mutations at positions 26 and 31 
(L26F/S31N), 27 and 31 (V27A/S31N) were also found in these hosts (Table 4). However, No 
G34E resistant mutation at position 34 was detected in avian influenza viruses. Similarly, Avian 
and swine influenza viruses also had no L38F adamantane resistant mutation at position 38. No- 
tably, adamantane resistance mutations in all hosts were detected most frequently at position 31 
(S31N), followed by position 27 (V27A). The most common drug resistance mutation was S31N 
in human, avian, and swine influenza viruses with the frequencies of 98.7% (10813/10954), 
88.1% (909/1032), and 77.4% (1595/2062), respectively. In addition, a small number of S31N 
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Table 4. Host distribution of adamantane-resistance mutations analyzed in this study. 


Adamantane-resistant 


The number of adamantane resistant mutations from all of the HA subtypes among different hosts 


Avian 


mutation 
Human Swine Wild Canine Cheetah Ferret Mink Cat equine bat environment unknown 
boar 
L26F U 2 
V27A 63 5 1 1 1 
A30T 2 1 
S31N 10813 1595 4 2 4 2 4 1 S S 2 
G34E 2 
L38F 1 
L26F /S31N 10 29 1 
V27A/S31N 50 429 1 1 4 1 
A30T/S31N 2 
A30V/S31N 2 
G34E/S31N 1 1 
L38F/S31N 1 
Overall 10954 


2062 «1 4 2 4 3 5 1 3 57 4 


doi:10.1371/journal.pone.0119115.t004 


1032 


adamantane resistance mutations were also detected in canine, cheetah, ferret, mink, cat, equine, 
bat, and environment influenza viruses (Table 4). In S1-S6 Tables we show more adamantane 
resistant variants carrying the S31N resistance mutation at position 31 in H3N2 and HIN1 
human influenza viruses after 2003. Adamantane resistant strains of avian influenza H5N1 vi- 
ruses have also been found more frequently since 2000. Our study further suggested that the fre- 
quency of emergence of adamantane-resistant strains varied among different hosts and 
adamantane resistance mutations appeared across different hosts. 


Geographical distribution of adamantane-resistant influenza variants 


To further understand the geographical distribution of adamantane resistance, the M2 protein 
sequences related to drug-resistant phenotypes were analyzed focusing solely on geographic lo- 
cation of influenza variants with adamantane-resistant mutations. The geographic analysis 
summarized in Table 5 indicated a substantially rising percentage of adamantane-resistant in- 
fluenza viruses isolated from the world including specific countries in Africa, Asia, Europe, 
America, and Oceania. Adamantane resistant influenza A H1 variants were distributed widely 
in the world. Of these variants, 52.1% (4572/8777) distributed in several countries in America 
(United States, Canada, Mexico, and Nicaragua) and 26.4% (2315/8777) were detected in Asia 
(China, Singapore, Thailand, Hong Kong, South Korea, Cambodia, Viet Nam, Japan, and Tai- 
wan). In addition, about one-fifth of the influenza H1 variants distributed in Europe (United 
Kingdom, Spain, Greece, Germany, Finland, Russia, and Netherlands) (Table 5). Adamantane 
resistant influenza A H3 variants were obtained mainly from 26 countries, showing a wide dis- 
tribution across the world. Of those, 34.8% (1427/4107) isolates from Asia and 56.4% (2316/ 
4107) variants from America were adamantane resistant. All influenza A H6 and H9 variants 
with adamantane resistance only distributed in Asia (China, Taiwan, and Hong Kong). For in- 
fluenza A H5 adamantane resistance variants, a great number of them were distributed in Asia 
including China, Thailand, Cambodia, Viet Nam, and Indonesia, while a small number of 
them were found in Africa (Egypt). Importantly, the largest proportion of Asian drug-resistant 
avian influenza viruses of H5 and H9 subtypes occurred in China (Table 5). Influenza A H7 
adamantane resistance variants distributed primarily in China of Asia and United States of 
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Table 5. Geographical distribution of adamantane-resistant influenza variants. 


Geographiclocation 


Africa 

Kenya 

Egypt 

Other Africa? 
Total Africa 
Asia 

China 
Singapore 
Thailand 
Hong Kong 
South Korea 
Cambodia 
Viet Nam 
Japan 
Taiwan 
Indonesia 
Other Asia? 
Total Asia 
Europe 
United Kingdom 
Spain 

Greece 
Germany 
Finland 
Russia 
Netherlands 
Other Europe" 
Total Europe 
America 
United States 
Canada 
Mexico 
Nicaragua 
Other America? 
Total America 
Oceania 
Australia 
New Zealand 
Total Oceania 
Overall 


The number of adamantane-resistant variants from different HA subtypes 


Hi 


52 


13 
65/8777(0.7%) 


662 
539 
314 
216 
64 
39 
70 
102 
135 


174 
2315/8777(26.4%) 


422 
389 

122 

101 

128 

105 

45 

355 
1667/8777(19.0%) 


3644 
335 
174 
141 
278 
4572/8777(52.1%) 


150 
8 
158/8777(1.8%) 
8777 


H3 


42 


61 
103/4107(2.5%) 


532 
166 
61 
255 
2m 
30 
91 
51 
129 


85 
1427/4107(34.8%) 


10 
15 


48 
2 


15 
99 
189/4107(4.6%) 


2122 
11 

35 

111 

37 
2316/4107(56.4%) 


40 
32 
72/4107(1.8%) 
4107 


a Ethiopia, Mali,Nigeria, Senegal, South Africa, Uganda. 
P Bangladesh, India, Iran, Israel, Kazakhstan, Kuwait, Kyrgyzstan, Laos, Macau, Malaysia, Mongolia, Myanmar, Pakistan, Philippines, Qatar, Saudi 
Arabia, Sri Lanka, Turkey, Turkmenistan, United Arab Emirates. 
° Belarus, Belgium, Czech Republic, Denmark, Estonia, France, Hungary, Ireland, Italy, Norway, Poland, Portugal, Serbia, Sweden, Switzerland. 
g Argentina, Bolivia, Brazil, Chile, Costa Rica, Cuba, Dominican Republic, Ecuador, El Salvador, Guam, Guatemala, Haiti, Peru, Puerto Rico. 


doi:10.1371/journal.pone.0119115.t005 
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H5 H6 


52 
3 
55/804(6.8%) 


119 4 


108 
22 


731/804(90.9%) 25/25(100%) 


1 
1/804(0.1%) 


13 

0 

3 
17/804(2.1%) 


804 25 


H7 


93 


2 
97/144(67.4%) 


2 
2/144(1.4%) 


45 


45/144(31.3%) 


144 


H9 


124 


68 
265/265(100%) 


265 
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Table 6. Geographical distribution of adamantane-resistant mutations analyzed. 


Geographiclocation The number of adamantane-resistant mutations from all of the HA subtypes 
L26F V27A A30T S31N G34E L38F L26F/ V27A/ A30T/ A30V/ G34E/ L38F/ 
S31N S31N S31N S31N S31N S31N 
Africa 
Kenya 94 
Egypt 52 
Other Africa“ Te 1 
Total Africa 223 1 
Asia 
China 1 35 8 1483 2 8 
Singapore 705 
Thailand 2 531 6 
Hong Kong 4 5 1 529 
South Korea 2 1 76 1 15 
Viet Nam 1 414 8 
Japan 1 151 ij 2 
Taiwan 285 1 1 d 
Indonesia 70 28 10 
Other AsiaP 5 16 437 3 5 
Total Asia 10 129 13 4639 1 13 49 1 
Europe 
United Kingdom 1 423 6 1 1 
Spain 4 360 40 
Greece 122 
Germany 137 12 
Finland 128 2 
Russia 152 1 
Netherlands 60 
Other Europe" 379 9 21 
Total Europe 5 1761 16 76 1 
America 
United States 1 5 1 5401 2 Ter 378 2 
Canada 1 342 1 1 
Mexico 1 221 
Nicaragua 251 1 
Other America® 330 4 3 
Total America 3 5 1 6545 2 21 382 1 2 1 
Oceania 
Australia 186 4 1 
New Zealand 3 1 36 
Total Oceania 3 1 222 4 
Overall 16 139 15 13390 2 1 51 511 2 2 2 1 


doi:10.1371/journal.pone.0119115.t006 


America, with the highest frequency of resistance found in China (64.6%). Similarly, geograph- 
ical analyses of adamantane-resistant mutations were also conducted in this study. The num- 
bers of adamantane-resistant mutations from all of the HA subtypes are shown in Table 6. The 
L26F mutation distributed primarily in Hong Kong and New Zealand, the V27A mutation 
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mainly distributed in China and Indonesia, and the A30T mutation distributed primarily in 
China. The S31N mutation distributed across the world (Table 6). This is the most common 
mutation known to confer resistance to adamantanes. In contrast, the G34E and L38F muta- 
tions only distributed in United States and South Korea, respectively. For double adamantane 
resistance mutations, the L26F/S31N mutations distributed mainly in Thailand, United King- 
dom, and United States. The majority of the V27A/S31N mutations emerged in United States, 
Spain, South Korea, and Indonesia (‘Table 6). The A30T/S31N mutations were found in Nicara- 
gua and Australia, and the G34E/S31N mutations occured in United Kingdom and Canada. 
The L38F/S31N and A30V/S31N mutations distributed only in Taiwan and United States, re- 
spectively. Our study further suggested adamantane resistance mutations appeared across dif- 
ferent geographic regions and the frequency of emergence of adamantane-resistant strains 
varied among different geographical areas. 


Temporal distribution of adamantane resistant influenza variants 


To better understand time distribution of adamantane resistant influenza variants, the M2 pro- 
tein sequence analyses were conducted focusing attention on years of isolation of influenza vi- 
ruses with adamantane resistance. The numbers and frequencies of adamantane resistant 
variants from H1N1, H3N2, H5N1, H6, H7, and H9N2 subtype influenza viruses among differ- 
ent years of isolation are given in Table 7-12. The trend for overall data showed a considerable 
increase. And the frequency of emergence of adamantane-resistant variants from different HA 
subtypes varied among different years of isolation during 2001-2013. For H1 subtype, tempo- 
ral distribution of adamantane-resistant influenza H1N1 variants during 1933-2013 revealed 
substantial increases in rates for USA, China, Singapore, United Kingdom, Spain, Canada, 
Thailand, Hong Kong, Mexico, Australia, Nicaragua, Finland, Greece, Russia and Taiwan 
(Table 7). Trends of rate changes from these countries were not considerable before the spikes 
occurred in 2009, which indicate the existence of such spikes and their considerable effects on 
the overall trends. As shown in Table 7, the highest level of adamantane resistance was contin- 
ually seen during 2009-2013 at 100% in many countries (e.g. Singapore, Spain, Portugal, Nica- 
ragua, Chile). For H3 subtype, temporal distribution of adamantane-resistant variants from 
H3N2 influenza viruses during 1968-2013 showed substantial increases in rates from many 
countries and regions across the five continents: USA, China, Hong Kong, Singapore, Taiwan, 
Nicaragua, Vietnam, Thailand, Uganda, Japan, Russia, Germany, Kenya, and Australia 

(Table 8). The frequency of adamantane-resistant H3N2 influenza viruses has increased sharp- 
ly since 2003, particularly in China. The proportion of isolates resistant to adamantanes contin- 
ued to increase during 2005-2006. Specifically, isolates collected during this period in 
Vietnam, Spain, Mexico, Cuba, and New Zealand exhibited 100% drug resistance. During 
2003-2009, resistants frequency increased considerably in mainland China, Hong Kong, Tai- 
wan and Thailand. During 2005-2013, there were substantial increases in the rates in USA and 
Singapore. In contrast, no adamantane resistance was observed in United Kingdom during this 
period. During 2008-2010, the overall rate change in Kenya, Uganda, Russia, Netherlands, Me- 
xico, and Japan was considerable. 

For the H5 subtype, temporal distribution of adamantane-resistant influenza H5N1 variants 
during 1959-2013 indicated considerable rate changes from six distinctly important countries 
including China, Thailand, Cambodia, Vietnam, and Indonesia of Asia and Egypt of Africa. 
Furthermore, rates among these countries were also considerably different (Table 9). For in- 
stance, the frequency of H5N1 adamantane-resistant viruses obtained from Indonesia began 
increasing in 2003, with a considerable spike (100%) between 2007 and 2008. This trend of 
having spikes was followed by an increasing frequency in Thailand and Vietnam in 2003. 
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Frequency and Distribution of M2 Gene Mutations 


Table 9. Temporal distribution of adamantane-resistant influenza variants from H5N1 subtype viruses during 1959—2013. 


Geographiclocation 


Egypt 

China 
Thailand 
Hong Kong 
South Korea 
Cambodia 
Viet Nam 
Indonesia 
Laos 

India 
Bangladesh 
Saudi Arabia 
Iran 
Malaysia 
Israel 
Elsalvador 


Belgium 


The number and frequency of adamantane-resistant variants from H5N1 subtype influenza viruses 


By 
2000 
0/10 
0/10 
0/0 
0/44 
0⁄1 
0/0 
0/74 
0/56 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 


0/0 


2001 


0/0 
2/18 


(11.1%) 
1⁄1 


(100%) 
0/6 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
0/0 
val 


(100%) 
0/0 


doi:10.1371/journal.pone.0119115.t009 


2002 


0/0 
7/32 
(21.9%) 
0/0 
2/15 
(13.3%) 
0/0 

0/0 

0/0 

0/0 

0/0 

0/0 

0/0 

0/0 

0/0 

0/0 

0/0 


0/0 


0/0 


2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 


0/0 0/0 0/0 0⁄4 2/16 18/23 13/16 15/17 1/10 0/2 0/0 
(12.5%) (78.3%) (81.3%) (88.2%) (10.0%) 
15/55 8/55 12/146 28/119 7/32 11/20 11/25 0/3 6/21 1⁄1 0/0 
(27.3%) (14.5%) (8.2%) (23.5%) (21.9%) (55.0%) (44.0%) (28.6%) (100%) 
3/3 85/85 37/37 WD 6/7 18/18 0/0 3/3 0/0 0/0 0/0 
(100%) (100%) (100%) (100%) (85.7%) (100%) (100%) 
2/7 0⁄1 0/0 0/18 1/27 2/6 0/0 0/2 0/3 0/0 0/0 
(28.6%) (3.7%) (33.3%) 
0/0 0⁄1 0/0 oe 0/0 0/6 0⁄4 07 1/25 0/0 0/0 
(4.0%) 
0/0 1/1 11/11 13/13 3/3 4/4 2/2 Gu 14/14 1⁄1 0/0 
(100%) (100%) (100%) (100%) (100%) (100%) (100%) (100%) (100%) 
17/17 8484 61/81 6/6 59/65 9/17 4/11 0/2 2/5 719 0/12 
(100%) (100%) (75.3%) (100%) (90.8%) (52.9%) (36.4%) (40.0%) (36.8%) 
5/9 2/9 17/24 58/76 21/23 5/5 0/0 0/0 0/0 0/0 0/0 
(55.6%) (22.2%) (70.8%) (76.3%) (91.3%) (100%) 
0/0 1⁄1 0/0 0/6 1/22 5/6 0/0 0/20 0/0 0/0 0/0 
(100%) (4.5%) (83.3%) 
0/0 0/0 0/0 3/8 on 0/30 0/10 3/3 0/8 0/0 0/0 
(87.5%) (100%) 
0/0 0/0 0/0 0/0 on 0⁄1 0/0 1/5 1/4 2/17 0/0 
(20.0%) (25.0%) (11.8%) 
0/0 0/0 Ko 0/0 2/3 0/0 0/0 0/0 0/0 0/0 0/0 
(66.7%) 
0/0 0/0 0/0 0⁄4 0/0 1⁄1 0/0 0/0 0/0 0/0 0/0 
(100%) 
0/0 0/0 0/0 0/0 1/1 0/0 0/0 0/0 0/0 0/0 0/0 
(100%) 
0/0 0/0 0/0 0⁄4 0/0 1⁄1 0/0 0/0 0/0 0/0 0/0 
(100%) 
0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 
0/0 1⁄1 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 0/0 
(100%) 


Notably, the percentage of adamantane resistant H5N1 viruses gathered from Cambodia is par- 
ticularly high and consistently reached 100% during 2004-2012. For H6 subtype, temporal dis- 
tribution of adamantane-resistant H6 influenza variants during 1999-2013 suggested that 
H6N1 and H6N2 mutants with considerable rate changes emerged mainly in Taiwan and Viet- 
nam, respectively (as shown in Table 10). No H6N2 resistant viruses were detected in samples 
collected from China. For H7 subtype, temporal distribution of adamantane-resistant H7 influ- 
enza variants during 1902-2013 suggested that H7N9, H7N7, and H7N6 mutants with consid- 
erable rate changes emerged mainly in China in 2013, whereas H7N3, H7N2, and H7N1 
mutants were respectively detected in Pakistan, USA, and Italy by 2005 (as shown in Table 11). 
For H9 subtype, temporal distribution of H9N2 adamantane resistant viruses during 1998- 
2013 displayed considerable changes in rates in Asia for specific countries or areas such as 
China, Hong Kong, Bangladesh, and United Arab (Table 12). In 2009, one or more drug-resis- 
tant virus was identified in five countries: China (20/23, 87.0%), South Korea (1/24, 4.2%), 
Vietnam (3/4, 75.0%), Bangladesh (4/7, 57.1%), and Iran (2/2, 100%). However, no resistance 
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Table 11. Temporal distribution of adamantane-resistant influenza variants from H7 subtypes viruses during 1902-2013. 


Geographiclocation 


The number and frequency of adamantane-resistant variants from H7 subtypes influenza viruses 


By 2000 
China (H7N9) 0/0 
Taiwan (H7N9) 0/0 
China (H7N7) 0/0 
Italy (H7N7) 1/3 

(33.3%) 
China (H7N6) 0/0 
Pakistan (H7N3) 0/15 
USA (H7N2) 11/99 

(11.1%) 
Italy (H7N1) 0/53 


doi:10.1371/journal.pone.0119115.t011 


2001 
0/0 


0/0 
0/0 
0/0 
0/0 
0/1 
15/32 


(46.9%) 
0/2 


2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 


0/0 0/0 0/0 0/0 on op Of oui op op op 63/63 
(100%) 

0/0 0/0 0/0 0/0 oo Of op op op op op 2/2 
(100%) 

0/0 07 0/0 0/0 on Off ou 03 os op 0 30/32 
(93.8%) 

0/0 0/0 0/0 on 02 o0 op op op op op op 

0/0 og 0/0 0/0 oo Of op 02 op op OO 1la 
(100%) 

0/0 0/2 2/4 0/0 oo Of op op op op op Of 

(50.0%) 

1/17 5/17 0/3 13/161 0/33 ou on op op op op 0/0 

(5.9%) (29.4%) (8.1%) 

1/1 0/0 0/0 0/0 oo Of op op op op op og 

(100%) 


was found in Hong Kong, Japan, Israel, and United Arab in the same period. By comparison, 
the frequency of resistant H9N2 variants isolated from China was consistently near 10% during 
2002-2005. 

In conclusion, our study further suggested the frequency of emergence of adamantane-resis- 
tant influenza strains varied among different years of isolation. 


Discussion 


Anti-influenza drugs play an important role in a comprehensive approach on controlling influ- 
enza A virus infections. Adamantanes had been excellent anti-influenza medicines until the re- 
cent emergence of resistant viruses. In this study, 31251 different subtype influenza A viruses 
(H1-H17) isolated in the world from 1902 to 2013 were assessed for resistance to adamantanes. 
On the basis of M2 protein sequence analysis, our study reveals a continuing worldwide in- 
crease in adamantane resistance, and suggests that the frequency of emergence of drug-resis- 
tant influenza variants varied among different HA subtypes, host species, years of isolation, 
and geographical areas. This study, which evaluates the frequency and distribution of 14132 
adamantane-resistant influenza viruses obtained worldwide, is the largest and most compre- 
hensive report on adamantane resistance to date. 

The frequency of resistance to adamantanes among circulating influenza A viruses has dra- 
matically increased over the past few years [16]. In this study we confirm that H1, H3, H5, H7, 
H9, and H17 subtype influenza A viruses with the specific resistance-associated mutations in 
their M2 genes exhibited high-level resistance to adamantanes. In contrast, the appearance of 
H2, H4, H6, H10, and H11 adamantane-resistant mutants was a rare event. However, no ada- 
mantane resistance viruses were identified among other HA subtypes (H8, H12-H16) 

(Table 1). Our findings extend those reported in 2005 [16] and 2006 [19] and document a con- 
tinuous and proportional increase in influenza A (H1, H3, H5, H7, H9) variants showing resis- 
tance to adamantanes. The global incidence of adamantane resistance among A (H3N2) has 
increased dramatically. This escalating trend in circulating H3N2 drug-resistant variants was 
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first observed among viruses isolated from Asia in 2000 and then in other regions of the world 
in 2005 (Table 8). A significant increase in resistance was also detected in many countries 
throughout the world among A (H1N1) viruses collected during 2009-2012 (Table 7). Al- 
though adamantane resistant H5N1 variants are present in Asia, their distribution appeared to 
be largely limited to Thailand, Cambodia, and Vietnam, It is worth noting that most H5N1 vi- 
ruses from Indonesia and China are sensitive to adamantine [5] (Table 9). The apparent geo- 
graphical disparity in the susceptibility of H5N1 isolates to adamantane is unexplained. 
Similarly, H9N2 viruses isolated from Asian countries including China, Hong Kong, Bangla- 
desh, and United Arab show significant increases in drug-resistance frequencies (Table 12). 
Particularlly, H7N9 resistant mutants only distributed in mainland China and Taiwan of Asia 
also show a significant increase in the incidence (100%) of adamantane resistance in 2003 
(Table 11). Our findings that epidemic and pandemic strains of influenza were identified first 
in Asia accord with the results from other previous studies [26]. 

Adamantane resistance in the influenza A virus is associated with six amino acid substitu- 
tions in the M2 protein according to previous reports [16]. The known mutations that confer 
adamantane resistance are L26F, V27A, A30T (A30V), S31N, G34E, and L38F [12, 22]. It was 
reported that mutated viruses may either lose the ability to bind M2 ion channel blockers, as 
with the S31N or A30T amino acid substitutions [27], or bind the blockers but retain M2 func- 
tion, as with amino acid replacements L26F or V27A at residue 26 or 27 [27]. In our study de- 
scribed here, the most common adamantane resistance mutation was S31N in influenza A 
viruses from human, avian, and swine (Table 4). Most adamantane-resistant influenza variants 
(95%) bear SN amino acid mutations at position 31 of the M2 protein, whereas only about 
1% have V—A mutations at aa 27. Drug-resistance mutations at other amino acid positions 
(L26F, A30T, G34E, and L38F) are extremely rare (their prevalence appeared to be < 0.2%) 
(Table 2). Among these mutations, the S31N is the most frequently reported [14, 19, 28], sug- 
gesting that variants containing the $31N substitution might possess a significant advantage on 
viral replication or transmission, leading to more efficient circulation. Additionally, the high 
levels of the dual L26I and S31N resistance mutations mainly detected in H5N1 variants ob- 
tained from Thailand, Vietnam, and Cambodia indicate that viruses carrying this dual motif 
are stably selected (S3 Table). Given that the dual L26I and S31N motif are identified firstly in 
China in 2002 in our study, it appears that China might have been the location of introduction 
or generation of viruses. Our results further suggest that M2 drug-resistant mutations could 
have occurred spontaneously before these drugs were developed, implying that greater caution 
is needed in the use of adamantanes. 

It is known that influenza viruses resistant to amantadine and rimantadine can emerge 
quickly when these drugs are being used to control influenza outbreaks [11, 29]. Adamantanes 
have been considered first-line drugs for the prophylaxis and treatment of influenza A virus in- 
fections. The widespread use of adamantanes has been associated with the rapid emergence of 
resistant viruses which are as genetically stable, virulent and transmissible as the wild-type 
virus [30]. Here we report the results of a comprehensive study focusing on adamantane resis- 
tance in influenza A viruses circulating worldwide. Our results reveal a significant increase in 
adamantane-resistance frequencies in influenza A H1N1, H3N2, H5N1, H6, H7, and H9N2 vi- 
ruses circulating during 2001-2013 (Table 7-12). The high levels of adamantane-resistance 
found in our study are comparable with other studies monitoring resistance worldwide [17, 18, 
19, 31]. For example, adamantane resistance in A (H3N2) strains continued to increase in 
many countries from 2005 to 2007 (e.g. New Zealand, USA, Malaysia, and China) [31] 

(Table 8). However, resistance in A (H1N1) strains had been quite variable during that time pe- 
riod in some countries such as Japan, Australia, USA and Canada [15, 32, 33] (Table 7). The 
high frequency of adamantane-resistant variants indicates that continuous global surveillance 
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and rapid identification of mutants are essential to monitor the emergence and spread of drug 
resistance, and to help with making informed decisions about antiviral usage in control of in- 
fluenza virus infections. 

'To our knowledge, this study is the first attempt to evaluate the frequency and distribution 
of adamantane resistant influenza variants throughout the world over a period of 111 years. Re- 
sults from our study provide further evidence that amino acid residues 26, 27, 30, 31, 34, and 
38 in the M2 protein play a major role in determining the resistance phenotype of adamantanes 
of influenza A viruses and highlight the necessity of monitoring the susceptibility of influenza 
A viruses to antiviral drugs. Our findings raise concerns about the increasing prevalence of ada- 
mantane-resistant influenza variants and draw attention to the importance of tracking the 
emergence and worldwide spread of drug-resistant variants. 
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